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Abstract. Cereal-legume intercrops represent a promising way of combining high productivity and agriculture sus- 
tainability. The benefits of cereal- legume mixtures are highly affected by species morphology and functioning, which 
determine the balance between competition and complementarity for resource acquisition. Studying species morpho- 
genesis, which controls plant architecture, is therefore of major interest. The morphogenesis of cultivated species has 
been mainly described in mono-specific growing conditions, although morphogenetic plasticity can occur in multi- 
specific stands. The aim of the present study was therefore to characterize the variability of the morphogenesis of 
pea plants grown either in pure stands or mixed with wheat. This was achieved through a field experiment that in- 
cluded three pea cultivars with contrasting earliness (hr and HR type) and branching patterns. Results show that 
most of the assessed parameters of pea morphogenesis (phenology, branching, final number of vegetative organs 
and their kinetics of appearance) were mainly dependent on the considered genotype, which highlights the import- 
ance of the choice of cultivars in intercropping systems. There was however a low variability of pea morphogenesis 
between sole and mixed stands except for plant height and branching of the long-cycle cultivar. The information pro- 
vided in the present study at stand and plant scale can be used to build up structural-functional models. These models 
can contribute to improving the understanding of the functioning of cereal-legume intercrops and also to the defin- 
ition of plant ideotypes adapted to the growth in intercrops. 
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Introduction 

In order to ensure agriculture sustainability, efforts have 
been made by researchers and farmers to reduce the use 
of fertilizers and pesticides. This challenges the mainten- 
ance of efficient and profitable agrosystems being able to 
face demographic growth. Because of their ability to fix 



atmospheric N2, legume species can improve the sustain- 
ability of cropping systems by helping to decrease the use 
of nitrogen fertilizers and favouring the diversification of 
crop rotations (Crews and Peoples 2004; Due et al. 2010). 
Also, seeds or forages of legumes are among the richest 
sources of proteins in crops, with a high nutritional value 
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for animals (Due et ai 2010). However, the potential prod- 
uctivity of legumes has not been reached, mainly because 
of a strong sensitivity of these species to lodging and foliar 
diseases (Ney and Corrouee 2005). This is in particular the 
case for pea (Pisum sativum), which is the main source of 
vegetable proteins in Europe. In this context, the increas- 
ing interest in growing cereal- legume intercrops (IC), 
such as wheat- pea mixtures, represents an alternative 
for reintroducing legume species in cropping systems. Sev- 
eral studies reported that these mixtures can provide high 
and stable yields compared with pure mono-specific 
stands (Ofori and Stern 1987; Jensen 1996; Corre-Hellou 
et ai 2006; Hauggaord-Nielsen et ai 2008). Such advan- 
tages result from a balance between complementary 
(e.g. separate root and canopy areas) and competition 
processes for light, water and nitrogen that occur between 
intercropped species. These complex interactions depend 
on the pedo-climotic conditions, agricultural practices and 
also on the morphology and functioning of the compo- 
nent species (Corre-Hellou et ai 2006; Lounay et ai 
2009; Louarn et ai 2010; Naudin et ai 2010; Barillot 2012). 

The latter point is mainly related to the choice of culti- 
vors, which therefore appears as a determinant factor of 
(i) the proportion of each component species at harvest 
and (ii) mixture productivity. Cultivors are usually discrimi- 
nated according to their earliness, sensitivity to diseases or 
potential yield. However, in the particular case of multi- 
specific stands, the above-ground architecture of a culti- 
var, given by its geometry, optical properties and topology 
of the phytoelements (Godin 2000), should also be token 
into account. Indeed, plant architecture defines the plant 
interface with biotic (e.g. with Mycosphaerella pinodes; 
Beasse et ai 2000; Le May et ai 2009) and abiotic factors 
(e.g. light; Ross 1981). In the cose of multi-specific stands, 
the complementarity between the architecture of the 
mixed species represents a crucial issue as it will deter- 
mine their respective ability to compete for light that in 
turn drives the production and allocation of biomass 
(Vorlet-Groncher et ai 1993; Sinoquet and Caldwell 1995). 

For pea, several genes involved in the development of 
the above-ground architecture hove been identified (for 
a review see Huyghe 1998). For instance, numerous ramo- 
sus mutants were described because of their altered 
branching behaviour (Arumingtyas et ai 1992). Plant 
height con also be altered through mutations mode on 
genes involved in internode growth (Kusnadi et ai 
1992). Genetic control of the compound leaf shape of 
pea has also been assessed and appears to be related 
to the UNIFOLIAJA gene (Gourlay et ai 2000). Precocity 
of pea cultivors has been shown to be regulated by 
genes (Hr and if) that control the sensitivity to photo- 
period for floral initiation and flowering (Murfet 1973, 
1975). These studies have promoted the breeding of 



several pea cultivors with contrasting architectures that 
therefore constitute as much as potential combinations 
for wheat-pea IC. Characterizing the morphogenesis (se- 
quence of developmental and growth processes leading 
to plant architecture) of these various pea genotypes is 
therefore of major interest for improving the manage- 
ment of intercropped stands. Several descriptors can 
be used to characterize pea architecture, the most com- 
monly used being those related to the leaf area and its 
spatial distribution as this strongly determines a plant's 
ability to compete for light interception. On a finer 
scale, both the amount and distribution of foliar area 
ore related to the number and geometry of stems and 
leaves produced during the initiation of each phytomer 
by the apex. A phytomer is defined as a basic unit re- 
peated along the stem and including an internode, a 
node, 0 leaf and one or several axillary buds (Gray 1849; 
White 1979). The sharing of resources within multi- 
specific stands also depends on the respective height 
reached by the component species (Sinoquet and Cald- 
well 1995; Schwinning and Weiner 1998; Louarn et ai 
2010; Barillot et ai 2011, 2012). Although the architec- 
tural parameters involved in the leaf area and height of 
plants are key factors of the mixture development, they 
hove been mainly described in mono-specific growing 
conditions (for a review see Munier-Jolain et ai 2005a). 
However, the morphogenesis of plants can be highly plas- 
tic when facing environmental variations; hence the 
question arises as to whether morphogenetic variations 
can occur between mono- and multi-specific stands. 

The aim of the present study was therefore to charac- 
terize the variability of pea morphogenesis grown either 
in pure stands or mixed with wheat. In order to have a 
large range of plant architectures and morphogenetic re- 
sponses, 0 field experiment was performed using three 
pea cultivors with contrasting growth habits. The growth 
and phenology of the pea cultivors were measured regu- 
larly during their growing cycle. This study provides infor- 
mation at both stand and plant scale in order to identify 
plant traits of interest that can contribute to the concep- 
tion of plant ideotypes. 

Methods 

Plant material and growing conditions 

A field experiment was carried out in 2010-11 at Broin- 
sur-l'Authion, western France (47°26'N, 00°26'W) in a 
clay soil (51 % clay, 26 % silt and 23 % sand). Daily mean 
air temperature, precipitation and photosyntheticolly 
active radiation (PAR) were recorded by a standard auto- 
matic agro-meteorological station located close to the field. 

Winter wheat (Triticum aestivum) cv. Cezanne and 
three cultivors of winter field pea (Pisum sativum), cv. 
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Lucy (hr type), AOPHIO (hr type) and 886/01 (HR type), 
were sown on 28 November 2010 in sole crops (SC) and 
IC. Unlike hr types, flowering of tlie HR cultivar is sensi- 
tive to pliotoperiod. Tlie sowing density of SC was 
250 plants m"^ for wheat. Optimal densities of pea culti- 
vors were chosen with respect to their ability for lateral 
development and the underlying risk of lodging. Sole 
crops composed of pea cultivors Lucy and AOPHIO were 
sown at 80 plants m"^ whereas cultivar 886/01 was 
sown at 40 plants m"^. Intercrops followed a substitutive 
design where the two species were mixed within the row. 
Wheat and pea crops in IC were sown at half their re- 
spective density in pure stands. From seedbed prepar- 
ation to harvest, local agronomic recommendations 
were followed and pest and weed were chemically con- 
trolled. Stands of sole wheat were fertilized with 14 g 
N m"^ whereas pea SC and wheat- pea IC were not sup- 
plied with external nitrogen. 

Statistical analyses described below were thus per- 
formed considering two factors: (i) pea genotype with 
three levels (Lucy, AOPHIO and 886/01) and (ii) cropping 
system with two levels (SC and IC). A sole crop of wheat 
was added to those six treatments but was not consid- 
ered for statistical analyses. These seven treatments 
(3 SC of pea, 3 IC and 1 SC of wheat) were arranged in ex- 
perimental units of 1.2 X 10 m^ within a randomized 
complete block design with three replicates. 



Wit and Von den Bergh (1965). Land equivalent ratio is 
the sum of partial LER values for wheat (LERwheat) and 
pea (LERpea): 



LERwheat — 7^ — . LER| 
'SCw 



pea - 



'ICp 
I^SCp' 



LER - LERwheat + LERpea 



where Vicw and Vscw ore yields of wheat and pea in IC, re- 
spectively, and Vscw and Vscp are yields of wheat and pea 
in SC, respectively. Land equivalent ratio values above 1 
indicate a benefit of intercropping over sole cropping. 

On the other hand, specific measurements on pea cul- 
tivors were made at plant scale. The morphogenesis of 
five pea plants per plot was characterized for each vege- 
tative axis, i.e. main stems and lateral branches. Only one 
branch at each nodal position of the main stem was fol- 
lowed up. Branches were denoted according to their 
topological position, i.e. main stems were denoted as 
Axis-0, then branches that emerged from node n of the 
main stem were referenced as Axis-n. For each axis 
group, the kinetics of phytomer appearance (unfolding 
leaf visible to the naked eye) were measured and fitted 
with Schnute's non-linear model (Schnute 1981) using 
the least-squares method. The model is written as 



-A(t) 



l/B 



Plant sampling and measurement of pea 
morphogenesis 

On the one hand, integrated parameters defined at can- 
opy scale (biomass, height, yield) were measured in each 
plot. Samplings were carried out on 0.75 m^ in the centre 
of each experimental unit. The above-ground biomass 
and the maximal height of each SC (wheat and pea) 
and IC plot were measured during the growing cycle at 
645, 1525 growing degree days (GDD) from emergence 
(base temperature 0 °C) and, lastly, at crop maturity 
(Table 1). The land equivalent ratio (LER) for grain yields 
of wheat-pea IC was also estimated according to De 



where y is the number of visible phytomers; estimated 
parameters ore A and 6 which implicitly define the 
shape of the curve; tmox is the last value of the time (t) do- 
main for which the model is fitted, corresponding to the 
end of the vegetative development of the stem; param- 
eter ymax is the value of Y at tmox and e is the residual. 
Parameters were optimized using the Levenberg- 
Marquardt iterative method (Marquardt 1963) with auto- 
matic computation of the analytical partial derivatives. 
The fitting procedure was performed for each axis group 
of each plant. The first derivatives of Schnute's adjust- 
ments were also used in order to estimate the rotes of 
phytomer production of the pea cultivors. 



Table 1. Harvest time, expressed In growing degree day (DD) from 
emergence (base, 0 °C), of pea and wheat grown In SC and In IC. 



Species 


Genotype 


Harvest time 
in SC (DD) 


Harvest time 
in IC (DD) 


Pea 


Lucy 


1900 


2275 


Pea 


AOPHIO 


1985 


2275 


Pea 


886/01 


2130 


2275 


Wheat 


Cezanne 


2275 


2275 



Statistical analyses 

Exploratory data analysis, analysis of variance and non- 
linear regression techniques were performed with R 
software (R Development Core Team 2012). Analyses of 
variance (ANOVAs) were performed following a two- 
factor linear model such that 

yijk = iJ-+Bi + Gj + Ck + {Gx Qji, + Eijk 

where y is any dependent variable, |jl the mean value ofY, 
6/ the effect of block /, Gj the effect of genotype j, the 
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effect of cropping system k (either SC or IC) and e the ran- 
dom error of measurement ijk. 

Normal distributions of tiie residuals of ANOVAs as well 
as those of Schnute's adjustments were tested using the 
Shapiro- Wilk test. Homoscedasticity was checked by ran- 
dom distribution of the residuals. Tukey's HSD tests were 
used for mean separation when three or more means 
were compared. 

Results 

Environmental conditions during crop growth are sum- 
marized in Fig. 1. The daily average air temperature ranged 
from -4 °Cto27 °Con31Januaryand27 June, respectively 
(Fig. lA). Irrigation supplied 30 mm of water on 21 April 
and 27 May 2011. The daily cumulated PAR (Fig. IB) 
ranged from 2.70 to 108.70 mol m"^ on 31 December 
2010 and 25 June 2011, respectively. 

Growth of SC and IC stands 

The results described in this section are derived from 
measurements at the whole stand scale. First, the bio- 
mass accumulation of each crop is shown in Fig. 2. The 



A 2000 




Sep Oct Nov Dec Jan Feb Mar Apr May Jun Jul 

Month 

Figure 1. Meteorological conditions during the growing season 
2010-11 at Brain-sur-l'Authion, France. Daily mean air tempera- 
tures and rainfall are shown in (A). Vertical arrows represent a 
water supply of 30 mm by irrigation. Daily cumulated PAR is shown 
in (B). The horizontal arrow represents the growing period of pea and 
wheat stands. 
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Thermal time ("C day) 



B 2000 




Figure 2. (A) Accumulation of above-ground biomass as a function 
of thermal time from emergence (base temperature = 0 C). Sole 
crops are shown in closed symbols and IC in open symbols. Pea cul- 
tivors are in solid lines: Lucy is denoted by circles (filled, open), 
AOPHIO by squares (filled, open) and 886/01 by diamonds (filled, 
open). Wheat is in dotted lines with triangles (filled) in SC, and 
with the corresponding pea cultivar symbol in IC (open circles with 
'Lucy', open squares with 'AOPHIO' and open diamonds with '886/ 
01'). (B) Contribution of wheat (white bars) and pea (block bars) in 
the final biomass reached by SC and IC stands. 



above-ground biomass of the two species increased 
from 645 to 1525 GDD and then slowed down until 
maturity (Fig. 2A). Wheat SC showed the highest amount 
of biomass during the growing cycle and finally reached 
1830 gm"^. The final above-ground biomass accu- 
mulated by the three pea cultivors in SC ranged from 
880 to 1275 g m"^. On average, '886/01' exhibited the 
lowest biomass in SC throughout the growing cycle and 
finally reached 880 + 67 g m"^. This lower biomass has to 
be related to its sowing density, which was 50 % of the 
other cultivors (see the Methods section). Intercropped 
wheat (IC stands) accumulated on overage 1230 g m"^ 
of biomass of all IC stands token together. Pea grown in 
IC stands produced 220-325 g m"^ of biomass at matur- 
ity. At this stage of development, the final biomass 
reached by IC stands was 1500 g m"^ averaged across 
pea cultivors (Fig. 2B). Wheat in IC contributed to the 
main port of the mixture biomass (on overage 3.75 
times the biomass of pea). Although cultivar 886/01 
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exhibited tlie lowest biomass in SC stands, this cultivar 
produced more biomass when intercropped with wheat 
compared with 'Lucy' and 'AOPHIO', despite being sown 
at half density. 

Harvest of IC stands (Table 1) was performed when 
wheat reached maturity (2275 GDD). For pea grown in 
SC, harvest was made earlier, i.e. between 1900 and 
2130 GDD depending on the cultivar. Cultivors Lucy and 
AOPHIO (hr types) reached their maturity earlier than 
'886/01'. Indeed, cultivar 886/01 (HR type) needs a longer 
photoperiod to reach flowering and was therefore har- 
vested later than hr cultivors. As a result, the maturity 
of cultivar 886/01 and that of wheat were reached in a 
similar period (2130 and 2275 GDD, respectively). 

Land equivalent ratios for groin yield (LERs) were esti- 
mated for IC (Fig. 3). Although wheat SC stands were fer- 
tilized with nitrogen and IC stands were not, the partial 
LER of wheat was systematically higher than 0.5 (0.65 
on average) whatever the companion pea cultivar. Thus, 
wheat yields in IC were higher than those in SC when nor- 
malized by the sowing density (wheat density in IC was 
half that of SC but yields were reduced by <35 %). Land 
equivalent ratios of IC stands composed of cultivors Lucy 
and AOPHIO were slightly < 1, meaning that the cumula- 
tive yield of wheat and pea in IC stands was lower than 
the sum of each SC yield. Analysis of partial LERs of 
'Lucy' and 'AOPHIO' (0.3 and 0.2, respectively) revealed 
that their yield markedly decreased in IC compared with 
pure stands. In contrast, the yield of IC based on cultivar 
886/01 was 25 % higher than the cumulative yield of 
wheat and '886/01' in SC. In these IC stands, the yield 
of '886/01' was strongly increased compared with SC con- 
ditions (LERsss/oi = 0.62 on overage). Interestingly, the 
ratio of the final above-ground biomass of '886/01' in IC 
to that of '886/01' in SC (0.37, Fig. 2B) was smaller than 



1.6 
1.4 
1.2 



<; 1.0 



a 0.8- 



2. 0.6 
0.4 
0.2- 
0.0 



i LER Wheal 
I LER Pea 




Lucy AOPHIO 886/01 

Figure 3. Land equivalent ratio of eacli wheat- pea mixture. Land 
equivalent ratio values >1 indicate a benefit of intercropping over 
sole cropping. N.B.: unlike pure stands of pea and IC, wheat SC 
were fertilized with external nitrogen. 



the groin yield ratio (LERsee/oi)- This means that the har- 
vest index (ratio groin yield/above-ground biomass) was 
particularly high in intercropped '886/01'. 

Wheat plants grown in SC and IC stands exhibited simi- 
lar height and reached a maximum of 0.95 m (Fig. 4). 
Wheat grown in IC was taller than pea from the early 
stages of development and whatever the pea cultivar. 
The height reached by the pea cultivors was variable, de- 
pending on the cropping system (SC or IC stand). The can- 
opy height of sole pea crops reached a maximum of 
0.83-0.94 m for 'Lucy' and 'AOPHIO', respectively. The 
height of sole pea finally decreased dramatically at the 
end of the growing cycle as a result of plant lodging. How- 
ever, the height of pea cultivors grown in IC stands re- 
mained over 0.70 m, meaning that pea plants grown in 
IC were staked by wheat stems thus preventing pea lodg- 
ing. The species height ratio (peo/wheot) in IC is shown in 
Table 2. The height ratio ranged from 0.33 to 0.48 at 700 
GDD, meaning that pea cultivors were strongly domi- 
nated in the first stages of development, especially for 




1000 1500 
Thermal time (°C day) 



2500 



Figure 4. Observed height of canopies during the growing cycle. Sole 
crops are shown in closed symbols and IC in open symbols. Pea cul- 
tivors are in solid lines: Lucy is denoted by circles (filled, open), 
AOPHIO by squares (filled, open) and 886/01 by diamonds (filled, 
open). Wheat is in dotted lines with triangles (filled) in SC, and 
with the corresponding pea cultivar symbol in IC (open circles with 
'Lucy', open squares with 'AOPHIO' and open diamonds with '886/ 
01'). 



Table 2. Height ratio of wheat-pea mixtures according to the 
growing degree day (GDD). 



Genotype 


Height ratio (pea/wheat) 






700 GDD 


1330 GDD 


1525 GDD 


2275 GDD 


Lucy 


0.45 


0.87 


0.85 


0.80 


AOPHIO 


0.48 


0.93 


0.93 


0.96 


886/01 


0.33 


0.97 


1.02 


1.03 
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'886/01'. From 1330 GDD on, the Lucy and AOPHIO culti- 
vors were slightly shorter than wheat, whereas '886/01' 
reoched the height of wheat. 

Morphogenesis of pea cultivors 

In a second step, measurements were made at plant 
scale in order to compare the morphogenesis of pea cul- 
tivors grown in pure stands with those intercropped with 
wheat. 

Lateral branching of pea plants. The total number of 
branches (including non-flowering ones) produced by 
individual pea plants is shown in Fig. 5A according to 
their nodal position on main stems. Most lateral 
branches emerged from the first and second phytomer 
of the main stems (Axis-1 and -2) whatever the cropping 
system (SC and IC) and pea cultivor. Very few branches 
were produced on the third phytomer of the main stems 
(only nine branches, all cultivors and cropping systems 
taken together). Cultivors Lucy and AOPHIO grown in SC 
developed on average 3-4 branches per plant whereas 
cultivor 886/01 in SC was the most branching with 
about 6 branches per plant on overage. Nevertheless, 
0 significant effect of the genotype on the number of 
branches was only found for Axis-2 (ANOVA F2,82 = 



A 5 




LucySC LucylC AOPHIO SC AOPHIO IC 886n]| SC 886/0110 

Figure 5. Lateral branching of pea cultivors (n = 15 plants for each 
condition). (A) Number of lateral branches developed by the pea cul- 
tivors Lucy, AOPHIO and 886/01 grown in SC and in IC. (B) Time of 
branching expressed in thermal time from crop emergence. 
Branches were distinguished according to their nodal position on 
the main stem. Axis-1: branches emerged at the first node; Axis-2: 
second node; and Axis-3: third node. 



13.93, P< 0.001). Indeed, cultivar 886/01 produced 
significantly more branches of type Axis-2 than Lucy and 
AOPHIO cultivors did (HSD P < 0.001). The ability of the 
886/01 cultivar to develop numerous lateral branches 
oppeored to compensate for its lower sowing density as 
shown by its biomoss accumulation, which was similar 
to those of cultivors Lucy and AOPHIO (Fig. 2). Sig- 
nificant effects of the cropping system on the number of 
Axis-1 (Fi,82 = 14.33, P < 0.001) and Axis-2 (Fi,82 = 5.48, P 
< 0.05) were olso found. Indeed, pea cultivors grown in IC 
tended to develop fewer branches than in SC. The most 
drastic decrease was observed for cultivar 886/01, which 
produced 40% less branches in IC (HSD P< 0.01 and 
<0.05 for Axis-1 and -2, respectively). However, for 
'Lucy' (-11 %) and 'AOPHIO' (-22 %) intercropped with 
wheat, it was not possible to detect any significant 
differences in the number of branches between the two 
cropping systems, probably because of the high variability 
observed. 

Lateral branches emerged between 275 and 420 GDD, 
all cultivors and cropping systems token together 
(Fig. 5B). Although they were initiated later, branches 
that developed on the second phytomer of the main 
stem (Axis-2) appeared 70 GDD earlier than those located 
on the first node (Axis-1). Indeed, in most coses, the first 
phytomer (corrying the first vestigial leaf) is located a few 
millimetres under the ground. This may mechanically 
delay the emergence of branches, depending for instance 
on the sowing depth or soil structure, which also offect 
the quantity of light perceived by the axillary bud (Jeudy 
and Munier-Joloin 2005). Statistical analyses showed 
that there was a significont effect of the genotype on 
the time of branch emergence (F2,87 = 5.09, P< 0.01; 
'^2,166 = 3.20, P < 0.05 for Axis-1 and -2, respectively). 
Indeed, branches developed on cultivar Lucy emerged 
significantly earlier than those of '886/01' (HSD P < 0.01 
and <0.05 for Axis-1 and -2, respectively), while 
'AOPHIO' hod on intermediary behaviour. There were 
however no significant differences in the time of branch 
emergence between SC and IC stands, meaning that 
this parameter of pea morphogenesis was mainly de- 
pendent on the genotype. 

As shown above, pea genotypes produced several 
branches during the first 500 GDD. This led to complex 
plant architectures with a high potential number of 
stems per plant. For the soke of clarity, we only consider 
hereafter one branch at eoch node of the main stems. 
As 0 result, plant architectures were simplified to a 
main stem potentially bearing one branch on its first 
three nodes (see Fig. 5A). However, only a port of these 
stems actually grew and completed their development 
up to flowering. As shown in Fig. 6, very few main stems 



6 AoB PLANTS www.aobplants.oxfordjournals.org 



© The Authors 2014 



Barillot et at. — Morphogenesis of pea cultivors in sole crops and intercropped witli wlieat 




Stem 

Figure 6. Frequency of pea plants grown in SC and in IC whose mea- 
sured branches have reached flowering. Branches were distin- 
guished according to their nodal position on the main stem. 
Axis-1: branches emerged at the first node; Axis-2: second node; 
and Axis-3: third node (n = 15 plants for each condition, only one 
single branch was considered at each nodal position). 

of pea flowered. This was in particular the case of cultivars 
886/01 (0 % of flowering main stems) and Lucy (0-20 % 
of flowering main stems for IC and SC, respectively). Flow- 
ering of main stems was significantly dependent on 
the cultivar {Fj^io = 6.83, P < 0.05). Indeed, the propor- 
tion of 'AOPHIO' plants that hod flowering main stems 
(25-45 % in IC and SC, respectively) was statistically 
higher than that of cultivar 886/01 (HSD P < 0.05). Never- 
theless, most flowering stems of the three pea cultivars 
were Axis-1 and -2 branches. For cultivars Lucy and 
AOPHIO, there were 85 % more of Axis-1 branches that 
carried on growing until flowering compared with Axis-2 



(regardless of the cropping system). In contrast, the pro- 
portion of flowering Axis-1 of cultivar 886/01 was similar 
to Axis-2. Pea plants grown in SC and those intercropped 
with wheat exhibited similar proportions of flowering stems. 

Rate of phytomer appearance and final number of 
phytomers. Phytomer production by stem apices 
followed sigmoid-type dynamics as illustrated in Fig. 7. 
These dynamics were fitted with Schnute's function. In 
order to reduce the variability on Schnute's parameters 
(Table 3), stems that stop growing before flowering were 
not taken into account in the fitting procedure. The first 
three parameters of Schnute's function. A, B and tmax. 
respectively, ranged from 0.36 to 4.30 x 10"^, 0.13 to 
0.56 and 1164 to 1476 GDD averaged across cultivars 
and cropping systems (SC and IC). The shape parameters 
A and B of Axis-2 branches appeared to be significantly 
dependent on pea cultivar (F2,33 = 6.43, P < 0.01; F2,33 = 
12.23, P < 0.001 for parameters A and 6, respectively). 
Indeed, parameters A and 6 of Axis-2 branches produced 
by cultivar 886/01 were statistically higher than those 
of 'Lucy' and 'AOPHIO' whatever the considered cropping 
system (HSD P < 0.05 and <0.01 for A and B, respectively). 
We did not find any significant effect of the cropping 
system on parameters A and 6. The duration of phy- 
tomer production (parameter tmax. expressed from 
stem emergence) appeared to be similar among the 
different cultivars and cropping systems. Supplementary 
statistical analyses were also performed in order to 
compare the kinetics of phytomer production of the 
different stems. For cultivar Lucy, the kinetics of phytomer 
production were statistically different between main 
stems and branches (HSD P < 0.05 and <0.01 for A and 
6, respectively). These analyses showed that parameter 

30 T 1 




Thermal Time (°C day) 

Figure 7. Typical kinetics of phytomer appearance on a vegetative 
stem of the pea cultivars. Observed values are in closed symbols. 
Non-linear adjustment (solid line) was performed using Schnute's 
equation. 
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Table 3. Parameters {A x 10~^, 6 and tmax) of Sclnnute's adjustments made on tine l<inetics of pliytomer appearance for eacli pea cultivar (Lucy, 
AOPHIO, 866/01) grown in SC or in IC. Sclinute's adjustments were performed for each flowering stem (0: main stem; 1: branch developed on the 
first node of the main stem; 2: second node). Indicated values are the mean + SD (n = 15 plants for each cultivar and cropping system). 



Axis Parameter Stand 







Lucy SC 


Lucy IC 


AOPHIO SC 


AOPHIO IC 


886/01 SC 


886/01 IC 


0 


A (xlO-^) 


0.36 + 0.18 




2.50 + 1.54 


0.82 + 1.16 








B 


0.56 + 0.04 




0.27 + 0.10 


0.45 + 0.17 








^max 


1415 + 148 




1378 + 204 


1257 + 82 






1 


A ixlQ-^) 


2.36 + 1.49 


2.80 + 0.99 


4.30 + 4.74 


2.48 + 2.18 


1.88 + 0.95 


2.01 + 1.89 




B 


0.19 + 0.12 


0.13 + 0.11 


0.18 + 0.22 


0.24 + 0.21 


0.16 + 0.13 


0.25 + 0.15 




^max 


1294 + 190 


1288 + 83 


1164 + 336 


1288 + 234 


1279 + 135 


1336 + 199 


2 


A (xlO-^) 


2.36 + 0.86 


2.63 + 0.53 


2.21 + 1.38 


2.49 + 0.071 


1.86 + 1.46 


2.18 + 1.27 




B 


0.18 + 0.09 


0.14 + 0.04 


0.19 + 0.09 


0.17 + 0.08 


0.34 + 0.08 


0.25 + 0.11 




^max 


1362 + 111 


1352 + 22 


1292 + 224 


1378 + 50 


1476 + 52 


1365 + 242 



Table 4. Rate of phytomer production of each pea cultivar (Lucy, AOPHIO, 866/01) grown in SC or in IC. Maximum rate of phytomer production 
(Vmaxi phytomer degree-day~^) and time at which it was reached (tv„„,, DD) are shown (mean + SD). Computations were carried out for each 
flowered stem (0: main stem; 1: branch developed on the first node of the main stem; 2: second node) (n = 15 plants for each cultivar and 
cropping system). 



Axis Parameter Stand 







Lucy SC 




Lucy IC 




AOPHIO SC 


AOPHIO IC 


886/01 SC 


886/01 IC 


0 


^max 


0.023 + 


0.002 






0.024 + 0.002 


0.033 + 0.006 










1284 + 


171.3 






641 + 148 


1004 + 279 






1 


^max 


0.026 + 


0.005 


0.030 + 


0.004 


0.056 + 0.072 


0.026 + 0.006 


0.050 + 0.020 


0.030 + 0.004 






965 + 


225 


836 + 


125 


1035 + 247 


950 + 259 


1141 + 170 


1121 + 224 


2 


^max 


0.026 + 


0.004 


0.030 + 


0.002 


0.034 + 0.013 


0.027 + 0.003 


0.030 + 0.001 


0.030 + 0.001 






811 + 


69 


761 + 


59 


929 + 100 


753 + 47 


1226 + 207 


1025 + 93 





















6 was significantly different between Axis-1 and -2 
branches of cultivar 886/01 (HSD P< 0.05). 

The maximum rate of phytomer production (Vmax) ond 
time at which it was reached (tv^„J were estimated by 
computing the first derivative of Schnute's adjustments 
(Table 4). The maximum rate of phytomer production 
ranged from 0.023 for 'Lucy' SC to 0.056 phytomer 
degree-day"^ for 'AOPHIO' SC, which means that at max- 
imum activity one phytomer appeared each at 20-45 
GDD. Parameter Vmax of each axis group was found to 
be dependent neither on pea cultivar nor on cropping sys- 
tem. The time of maximum rate of phytomer appearance 
(^v^x) was reached between about 640 for 'AOPHIO' SC 
and 1285 GDD for 'Lucy' SC. Parameter tv^^, of flowering 
main stems was significantly higher for 'Lucy' than for 
cultivar AOPHIO (HSD P < 0.05). Pea cultivors were also 
statistically different for parameter tv™, of Axis-1 and -2 
branches (Fz.ss = 5.84, P < 0.01; F2,32 = 37.43, P < 0.001, 



respectively). Compared with 'Lucy', the maximal rate of 
phytomer appearance of cultivar 886/01 was indeed 
reached significantly later for Axis-1 and -2 (HSD P< 
0.01 and <0.001, respectively). Parameter tv^„, of cultivar 
886/01 appeared to be significantly higher than that of 
'AOPHIO' but only for Axis-2 branches (HSD P < 0.001). 
For this axis group, tv^^^ was also found to be higher in 
SC (1020 GDD) than in IC (845 GDD; HSD P < 0.05). 

The last parameter of Schnute's adjustments to be ana- 
lysed is the final number of phytomers (ymax)- The num- 
ber of phytomers produced on flowering stems ranged 
from 9 to 32 averaged across pea cultivors and cropping 
(Fig. 8). A marginal proportion of flowering stems exhib- 
ited <15 phytomers (1 % for 'AOPHIO', including main 
stems and Axis-2 branches). Stems with >15 phytomers 
were mainly branches developed on the first and second 
phytomer of main stems. This behaviour was however dif- 
ferent in the case of cultivar AOPHIO, for which ~20 % of 
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Figure 8. Frequency of flowering stems according to the final number of phytomers. Measurements were mode on three pea cultivors: Lucy, 
AOPHIO and 886/01 grown in SC and in IC with wheat. Results are shown for stems that reached flowering. Main stems ore denoted as Axis-0 and 
branches were distinguished according to their nodal position on the main stem. Axis-1: branches developed at the first node; Axis-2: second 
node; and Axis-3: third node (n = 15 plants for each condition). 



flowering stems were main stems with >20 phytomers. 
Statistical analyses also showed that the final number 
of phytomers measured on Axis-1 and -2 branches 
was dependent on the cultivar (F2,6o = 58.24, 
P < 0.001; F2,37 = 25.66, P < 0.001 for Axis-l' and -2, re- 
spectively). Indeed, cultivar 886/01 produced significantly 
more phytomers than 'Lucy' and 'AOPHIO' did (HSD 
P < 0.001 for both Axis-1 and -2 branches). The final num- 
ber of phytomers was not statistically different for pea 
plants grown in pure stands or intercropped with wheat. 



Flowering. The reproductive development of pea cultivors 
was characterized by the time of flowering of each 
stem, as well as the nodal position of the first flower. 
Flowering stage was reached between 950 and 1400 
GDD for AOPHIO IC and 886/01, respectively (Fig. 9). For 
a given cultivar and cropping system, flowering was 
synchronized between main stems and branches. 
However, flowering of Axis-1 and -2 branches was 
significantly different among the pea cultivors (F2,6i = 
1717.29, P < 0.001; Fz.ay = 1452.78, P < 0.001 for Axis-1 
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Figure 9. Time of flowering of pea cultivors Lucy, AOPHIO and 886/01 grown in SC or in IC. IMoin stems are denoted as Axis-0 and brandies were 
distinguislied according to tlieir nodal position on tlie main stem. Axis-1: branches developed at the first node; Axis-2: second node; and Axis-3: 
third node (n = 15 plants for each condition). 



and -2, respectively). Indeed, the flowering stage of 
cultivar 886/01 was reached significantly later (1400 
GDD) than for 'Lucy' and 'AOPHIO' (HSD P< 0.001 for 
both cultivors). The first flowering phytomer (Fig. 10) 
was located between the ninth and the 28th phytomer 
averaged across pea cultivors and cropping systems. 
First flowers of cultivar 886/01 (Axis-1 and -2 branches) 
were observed at higher phytomer positions (24th 
phytomer on average) than those measured for 'Lucy' 
and 'AOPHIO' (HSD P< 0.001 for both axis group and 
cultivors). As observed for the time of flowering, the 
nodal position of the first flower was similar between 
pea plants grown in pure stands and those mixed with 
wheat whatever the cultivors. 

Discussion 

The aim of the present study was to address the question 
of the morphological responses of pea to the competition 
when intercropped with wheat. To this end, a field experi- 
ment was conducted on three pea cultivors that were 
characterized at both stand and plant scale when 
grown in pure stands and mixed with wheat. 

Although the present study was conducted throughout 
1 year only, our results on crop biomass, species height 
and maturity were consistent with previous studies per- 
formed on different cultivors grown under contrasting 
pedo-climatic conditions (e.g. Corre-Hellou et al. 2009; 
Noudin et al. 2010). The results presented in this study 



can therefore be assumed to be representative of 
the conditions commonly encountered in wheat- pea 
mixtures. These results were obtained from unfertilized 
mixtures as usually performed in cereal-legume inter- 
cropping systems. Furthermore, some authors (Jensen 
1996; Corre-Hellou et al. 2006) found that the contribu- 
tion of the component species to the biomass of the mix- 
ture was dependent on the available nitrogen, an 
increase of which enhances the growth of the cereal spe- 
cies. The level and timing of nitrogen fertilization (Noudin 
2009; Noudin et al. 2010) therefore constitute o key factor 
enabling one to manage the hierarchy between the 
mixed species. The present study also shows that the 
three pea cultivors exhibited a similar level of biomass, 
especially when grown in mixtures. This suggests that 
despite the genotypic differences and contrasting initial 
sowing densities, the three pea cultivors hod a similar 
overall development even when they were in competition 
with wheat. However, this does not necessarily mean that 
the morphological processes of the pea cultivors had 
similar responses to the competition with wheat, but 
their integration at the stand scale leads to an equivalent 
growth in biomass. 

Our results also show that pea was strongly affected by 
lodging when grown in SC. Corre-Hellou et al. (2011) re- 
ported that the high sensitivity of pea to lodging caused 
significant yield losses as well as on enhanced growth of 
weeds. Pea lodging was however strongly decreased in 
mixed stands, pea branches being stacked by wheat. 
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Figure 10. Plant frequency according to the first flowering phytomer of stems. Results are shown for the pea cultivors Lucy, AOPHIO and 886/01 
grown in SC or in IC (n = 15 plants for each condition). Main stems are denoted as Axis-0 and branches were distinguished according to their 
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Intercropping cereals and legumes is therefore a promis- 
ing way of both reintroducing legume species within 
agrosystems and solving the problems encountered in 
pure stands of legumes. Moreover, the height reached 
by each species in the canopy is an important feature of 
the stand which determines, but also emerges from, the 
competition processes occurring between plants. Compo- 
nent species height ratio has been widely shown to affect 
light sharing in a mixture (Sinoquet and Caldwell 1995; 
Louarn et al. 2010; Barillot et al. 2011, 2012) and is there- 
fore a strong component of the inter-specific competition 
occurring within the mixture. The results described in the 



present study illustrate that the species height ratio is not 
constant throughout the growing cycle; pea cultivors 
were much shorter than wheat until 700 GDD and then 
reached a similar height. Differences among pea cultivors 
were also observed but were not constant over time. 
Although '886/01' was the shorter one in the early stages 
of development, this cultivar finally reached the same 
height as that of wheat afterwards. Therefore, it seems 
that the competition which occurs among the inter- 
cropped species cannot be assessed by punctual mea- 
surements of the species height ratio (in particular 
during the early stages of development). This is 
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consistent with the findings reported in a previous study 
(Barillot et al. 2012) where a virtual plant approach was 
used to demonstrate that the ability of plants to intercept 
light was mainly determined by the architectural para- 
meters involved in (i) the LAI (number of branches and 
phytomers, leaf area) during the early stages and (ii) 
plant height (internode length, number of phytomers) 
once canopy closure was established. 

The time lag between the physiological maturity of 
wheat and pea is a well-known issue of these mixtures. 
The choice of harvest timing is indeed complicated by 
the fact that pea generally reaches its maturity earlier 
than wheat. Nevertheless, physiological maturity of pea 
varies among the cultivors according to their eorliness, 
which is assumed to be mainly driven by the sensitivity 
of flowering to the photoperiod that involves the Hr 
gene (Murfet 1973). The maturity of the HR cultivar 
(886/01) was therefore almost synchronized with that of 
wheat, whereas the hr cultivors had to be harvested earl- 
ier. Gaps of maturity, as encountered with hr cultivors, 
represent a strong practical constraint at harvest (Louarn 
et al. 2010); HR pea cultivors therefore appeared to be 
well suited to intercropping with wheat. 

In order to deepen the analysis of the variability of 
pea morphogenesis in response to intercropping, a com- 
parison was performed at plant scale with particular at- 
tention to pea branching, flowering, final number of 
phytomers and their kinetics of appearance. Branching 
has been shown to be dependent on several factors 
such as genotype, hormonal balance, environmental fac- 
tors, e.g. low temperatures (Jeudy and Munier-Jolain 
2005), and also plant density (Spies et al. 2010). In the 
present study, contrasting abilities for branching were in- 
deed found between the genotypes [Lucy-AOPHIO] and 
886/01. Cultivar 886/01 was the most branching cultivar, 
which balanced its lower sowing density (50 % less than 
'Lucy' and 'AOPHIO'). Moreover, the number of branches 
tended to decrease in IC compared with pure stands, in 
particular for cultivar 886/01. Some authors like Casal 
et al. (1986), Ballore and Casal (2000) or Evers et al. 
(2011) showed that branching of several species is af- 
fected by the quantity (PAR) and quality (red/for-red 
ratio) of light perceived by the axillary buds. In the present 
study, we can therefore hypothesize that quantity of light 
and/or its quality were quite similar between the respect- 
ive pure stands and IC of cultivors Lucy and AOPHIO. This 
would mean that the replacement of a 'Lucy' or an 
'AOPHIO' plant by o wheat one leads to similar variations 
of light microclimate. This could be the result of small dif- 
ferences in the architectural patterns of the two species 
in terms of leaf area, height, geometry and/or optical 
properties. As cultivar 886/01 has a late development 
(HR type), we can also hypothesize that when 



branching started, wheat plants were more developed 
than neighbour '886/01' pea plants would have been in 
0 pure stand. This could cause variations of the microcli- 
mate perceived by axillary buds, leading to on inhibition 
of branching. 

The kinetics of phytomer appearance were assessed for 
main stems and a randomly selected branch at each 
node by using non-linear fittings. Our analysis showed 
that there were few statistical differences between 
the parameters belonging to the different genotypes 
and cropping systems. It was only found that (i) Axis-2 
branches of cultivar 886/01 had kinetics different from 
those of 'Lucy' and 'AOPHIO' and (ii) the maximum rate 
of phytomer appearance of '886/01' was reached later 
compared with the other cultivors. These results mean 
that the kinetics of phytomer production of different 
stems can be analysed/modelled by using similar 
Schnute's functions, at least for Lucy and AOPHIO culti- 
vors whether they were grown in sole stands or mixed 
with wheat. Turc and Lecoeur (1997) also reported similar 
rates of leaf primordium initiation and emergence for 
contrasting plant growth rates, cultivors and sowing 
densities in spring pea. One drawback of using Schnute's 
function lies in the fact that some of the parameters, es- 
pecially A and 6, cannot be directly related to a biological 
meaning. It would be tempting to use linear regressions 
because of the reduced number of parameters and easy 
interpretation. However, phytomer production is not in- 
trinsically constant and is actually characterized by o 
maximum rote (which can be estimated by the derivative 
of Schnute's functions) and o time at which development 
stops. These aspects cannot be handled by linear models. 
Supplementary statistical analyses (data not shown) 
showed that the residual sum of squares was significantly 
higher for linear regressions than that obtained with 
Schnute's function (HSD P < 0.001). These tests also indi- 
cated that the residuals of most linear regressions were 
not normally distributed and have means differing from 
zero. Nevertheless, the estimated parameters derived 
from Schnute's adjustments were highly variable. This 
variability is related to pea branching which is rather com- 
plex, particularly in the case of winter-sown cultivors. 
Winter conditions often cause frost damage, which in- 
duces the cessation of the development of the main 
stems and the initiation of numerous branches at differ- 
ent times (Fig. 5A and B). The result is a high variability in 
the characteristics of branches. 

In the present study, a significant difference was ob- 
served among the pea genotypes for the final number 
of phytomers reached on stems. Indeed, '886/01' (HR 
type) was found to produce more phytomers than the 
other cultivors. Similar results were also reported for 
this particular cultivar but grown in controlled conditions 



AoB PLANTS www.aobplants.oxfordjournals.org 



© The Authors 2014 



Barillot et al. — 



Morphogenesis of pea cultivors in sole crops and intercropped witli wlieat 



and individual pots (Barillot et al. 2012). In contrast, tlie 
number of initiated phytomers was similar whether pea 
plants were grown in pure stands or mixed with wheat, 
whatever the genotype. Moreover, our results show that 
the canopy of the three pea cultivors was mainly com- 
posed of branches as main stems had stopped growing 
with few phytomers. As reported by Jeudy and Munier- 
Jolain (2005), the development of branches is increased 
in winter pea cultivors because of the frost damage ex- 
perienced on the apex of the main stem. Such conditions 
were encountered during the first months of the growing 
cycle (December- February; Fig. 1) which corresponds to 
the emergence of the lateral branches (Fig. 5B). 

Flowering is a crucial stage of the growing cycle that 
has been widely studied and used in order to model pea 
growth. Truong and Duthion (1993) showed that the time 
of flowering is a function of leaf appearance rate and pos- 
ition of the node bearing the first flower. The reproductive 
development of pea cultivors was therefore characterized 
by two main indicators: the nodal position of the first 
flower and its emergence time. As also reported by Jeuf- 
froy and Sebillotte (1997), we found similar time of flow- 
ering between main stems and basal branches (although 
these were produced later) for all cultivors and cropping 
systems. Furthermore, the position of the first flowering 
node was similar among the genotypes and cropping sys- 
tems. Some authors also showed that for a given geno- 
type, the position of the first flowering node was 
constant over various conditions (Roche et al. 1998; 
Munier-Joloin et al. 2005b). 

Finally these results highlight that in the present ex- 
periment, the morphogenesis of pea was mainly deter- 
mined by the genotype and was only little affected by 
the competition with wheat. This suggests that the archi- 
tectures of pea and wheat may be quite similar, so that 
the environmental conditions perceived by plants in the 
canopy (phylloclimate; Chelle 2005) were not strongly dif- 
ferent between sole pea crops and wheat- pea mixtures. 
Functional-structural models (Vos et al. 2010; DeJong 
et al. 2011) ore able to take into account the explicit archi- 
tecture of plants and its interactions with physiological 
processes and environmental conditions. Such models 
therefore constitute suitable tools for assessing these hy- 
potheses and can in particular be used to characterize the 
microclimate perceived by plants located in mono- and 
multi-specific stands. 

Conclusion 

To our knowledge, the present study is the first to compare 
the morphogenesis of pea grown in sole stands with that 
of pea grown intercropped with wheat. On the one hand, 
the present results show that most of the assessed 



parameters of pea morphogenesis (phenology, branching, 
final number of phytomers and their kinetics of appear- 
ance) were mainly dependent on the considered genotype. 
This emphasizes the importance of the selection of culti- 
vors, in particular for intercropping systems, as this will de- 
termine the level of competition and complementarity 
between the component species. On the other hand, 
there was a low variability of pea morphogenesis between 
sole and mixed stands except for plant height and branch- 
ing of the late cultivor 886/01. Complementary studies on 
wheat-peo mixtures under contrasting levels of nitrogen 
fertilization are now needed to provide information on 
how nitrogen would affect plant morphogenesis and inter- 
specific competition. The information provided in the pre- 
sent study can be used for modelling pea morphogenesis 
in pure and mixed stands and therefore contributes to a 
better understanding of the functioning of cereal -legume 
IC. This kind of approach is also well suited for the identi- 
fication of plant traits to be integrated in the definition of 
plant ideotypes. 
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